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Abstract 

The states of hydrogen atom with principal quantum number n < 3 and zero magnetic 
quantum number in constant homogeneous magnetic field TL are considered. The perturbation 
theory series is summed with the help of Borel transformation and conformal mapping of the 
' (— i ', Borel variable. Convergence of approximate energy eigenvalues and their agreement with 

Q_i| corresponding existing results are observed for external fields up to n 3 TL ~ 5. The possibility 

of restoring the asymptotic behaviour of energy levels using perturbation theory coefficients 
is also discussed. 
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The own magnetic fields of some astrophysical objects reach very high values 0- If we 
J>v are interesting of the atomic spectra in these external fields, it is convenient to introduce natural 
measure of field strength - the atomic magnetic field = e 3 m 2 c/h 3 = 2.55 x 10 9 G. The fields 
H up to one half of Ho are detected in vicinity of some white dwarfs. Neutron stars possess fields 
| up to ~ 10 4 7io- For correct interpretation of the observations results it is desirable to know the 
atomic hydrogen spectrum in this range of external fields. For this aim, computations based on 
adiabatic approach with Landau level as initial approximation were accomplished |J. It will be 
shown here which part of the desired external field range could be covered with the help of the 
usual expansion in powers of 7i, starting from the Coulomb levels of hydrogen atom. We involved 
C\| \ in the computations many orders of perturbation theory (up to 75th order). Summation of the 
series was performed with the help of Borel transformation, supplemented by conformal mapping 
q ■ of Borel variable. 

The Borel summation method was introduced into quantum field theory long enough (see 
e.g. ||). It has been tested on some quantum- mechanical problems (one of many examples is 
described in ||) and continues to find applications in modern works ||. Large hopes on the possi- 
bility to advance into strong coupling region were related with Borel summation of the perturbation 
series. Some rather simple problems, in which details can be traced and compared with correspond- 
ing exact results, supported this optimism. For example, for the funnel potential, V(r) = —1/r+gr, 
by applying conformal mapping of the Borel variable and Pade-summation of the Borel transfor- 
mant, the ground state energy at g —>■ oo was obtained in the form E(g) = Cg v with ~ 0.2 % 
precision for index v and ~ 5 % precision for coefficient C 1. It became clear later that such a 
successful summation presents a special but not the general case. One can guess that this success 
is a consequence of simplicity of this problem. In contrast, asymptotic behaviour of energy levels 
in Stark and Zeeman effects comes into action at very large external fields values. For the Stark 
effect it is practically impossible now to reach the region of true asymptotic by perturbation series 
summation. An intermediate linear asymptotic is observed instead [[?]. § . 
To introduce notations and scale we write down the Hamiltonian 

H = -V - - + ^(r 2 - z 2 ) =H + gK (1) 
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Here g = 7i 2 , and hereafter we use units h = c = m = e = l. In (1) we drop the elementary 
contribution of electron's spin and consider only states with magnetic quantum number m = 0. 
We can expand E(g) as a formal series in powers of g: 



E(g) = J2E k g k - (2) 



k=0 

Now, we have to obtain hypersucceptibilities E^. We could use the moment method for this aim. 
This method is especially useful in the cases when variables in the Schrodinger equation can not be 
separated. Obviously the Zeeman effect presents just such a problem. In the previous work || the 
moment method was applied to recurrent evaluation of hypersucceptibilities. Somewhat different 
version of the moment method was introduced in the work JTOj. 

For the four lower "isolated" hydrogen levels we immediately use here the results of ref. 0. Un- 
fortunately the computer code, employed in the work || for the relatively more complicated case 
of degenerate 3s and 3d states, contained a mistake?]. Therefore we carried out new computation of 
3s and 3d hypersucceptibilities. Results of computation for some orders are presented in Table 1. 



These results are in agreement with results of the work [TI| where high-order hypersucceptibilities 
were obtained at the first time (but the method used in [JXTJ] is much more complicated than the 
moment method is). 



Table 1 



Hypersucceptibilities of degenerate states 



k Ek for 3s state 



1 

2 
3 
4 
5 

25 
50 
75 



19.57851476711195477229924488394 

-7992.558488642566993349104381687 

9951240.466276842310264046307800 

-20931559882.53444368634980579917 

58826900682409.79349115290157121 

1.3793233851820609414463787913215xl0 94 

-9.3227132696889616617788676903516xl0 211 

2.8053533970811704326574930831176X10 340 



k 



Ek for 3d state 



1 

2 
3 
4 
5 

25 
50 
75 



5.171485232888045227700755116050 

-1017.425886357433006650895618312 

738127.8247387826897359536921995 

-923576528.5544112941189442008231 

1677908319019.727217770438272530 

1. 0431217771 758614011812311858395xl0 92 

-6.0721978561446884300072726553011 xlO 209 

1.7302552995055432680731087635037xl0 338 



As the order k increases, hypersucceptibilities grow as a factorial 

E k ^E k = {-l) k+l C nl a k n T{2k + (3 nl ), 



(3) 



1 We are thankful to Prof. V. D. Ovsyannikov for drawing our attention to this mistake. 
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where a n = (n 2 /ir) 2 , [3 n i = 2n — l + ^^-, and C n i are not essential for us; one can find values of them 
in || and references therein. Eq. (3) implies that series (2) is asymptotical and the formal sum of 
such a series is ambiguous. But in fact the choice of the summation method is restricted: physical 
considerations impose analytical properties of the function E(g), which the true sum of series (2) 
is to reproduce. In the unphysical region, at g < 0, the diamagnetic perturbation gH\ changes its 
sign, the total Hamiltonian becomes "open" and the possibility of a spontaneous ionization of the 
atom appears. Therefore energy eigenvalue should have imaginary part at g < and the function 
E(g) should have a cut along real negative semi-axis of g plane. Summation with the help of Borel 
transformation results in a function having the left cut, besides the discontinuity on this cut is a 
smooth function of g. 

The Borel transformant B(w) of function E(g) is a series with coefficients B^ = Ek/T(2k + bo): 

oo 

B{w) = Y,B k w\ (4) 

fc=0 

where bo is an arbitrary constant. The choice of b can affect, in principle, on the numerical results, 
but changing of its value within interval ~ [0.5 < bo < 5] reveals weakly, so the choice of bo was 
made rather by convenience. The numerical calculations in this work were performed at b = 3. The 
series (4) converges, as usual, within the circle \w\ < l/a n . It is easy to check that the singularity 
of B{w) is located at w = —l/a n , substituting asymptotical coefficients Ek in place of E^. Energy 
of the level is related with the function B{w) by an integral transform 

poo 

E(g) = / e- x B(gx 2 )x b °- l dx. (5) 
Jo 

For the numerical integration in the right hand side to be successful, an analytical continuation of 
B(w) from its convergence circle on the domain, containing the image of the entire real positive w 
semi-axis is required. For this aim we performed conformal mapping of the Borel variable w. Many 
sufficiently effective versions of this mapping are appropriate. The main point is that the nearest 
singularity of the Borel transformant should be removed to infinity. Here we used the mapping 

a n w 

y — 

1 + a n w 

which was employed in the work 0. As is explained in ||, this transformation is optimal in the 
sense that it diminishes the influence of all possible singularities of B(w) from the unphysical region. 
Transformation (6) is equivalent to the following series rearrangement 

B W = £ 1W, D„ = B , D m = £ } " - >1 W J , m > 1. (7, 

m=0 k=l v ' v ' 

To improve the convergence we applied Pade summation to rearranged series (7) 

B{w) « [M/N](y) = P M {y)/Q N {y), (8) 

where Pm and Qn are polynomials of degree M and N respectively. 

We performed computations using various Pade approximants and straightforward summation 
of the rearranged series (7). To illustrate the influence of computational accuracy on summation 
results we compared ones made in double precision (16 decimal digits) with these in quadruple 
precision (32 decimal digits). 

Some graphs of the obtained binding energy £(TC) = ~Tt — E(Ti. 2 ) as a function of parameter 

7 ee n 3 7i are given in Figs. 1-3. As compared with the previous work 0, the region of external field 
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values for which these eigenvalues are successfully recovered is extended by a factor of about 5. As 
usual the precision of the sum considerably increases at lower Ti values. For instance in the case of 
3d state with the help of approximants [N/N](y) in the range 27 < N < 32 we get binding energies 
with 4 stable decimal digits at 7 = 4, with 6 digits at 7 = 3 and with 12 digits at 7 = 1. Note 
that in the work [[J Pade approximants were applied immediately to summation of divergent series 
(2). These approximants imitate the discontinuity on the cut g < by a set of delta-functions, 
and it is a very rough approximation. At the same time as a result of Borel summation the same 
discontinuity is represented by a smooth function of g. Our calculations confirmed that mapping (6) 
is indeed very effective: after this mapping Pade summation of the Borel transformant improves the 
convergence only a little and for some cases its straightforward summation appears to be sufficient 
- see Figs. 1-3. 

One technical detail is of principal importance for perturbation series summation by any method. 
The precision of the entire chain of computations must increase as the number of involved successive 
terms increases. This is simply a consequence of the fact that the sum, being of the order of unity, 
arises as a result of very large alternating sign terms compensation. 

It seems at first sight that the requirement of high precision is not necessary for the Borel 
transformant: all essential alternating sign coefficients Bk have about the same order. But any 
numerical procedure of analytical continuation usually requires high precision. Turning to series 
rearrangement (7) we see that binomial coefficients entering the sum for Dk are changing 20 orders of 
magnitude (in the present case). Obviously, we have enormous loss of precision performing the sum 
for Dk in (7). Therefore, if we want to use all Bk up to 75th order, the precision of Bk coefficients 
should be better than ~ 10 -20 . In our calculations, the precision of Ek and, consequently, the 
precision of Bk was ~ 10~ 30 , so the precision of Dk was decreasing from 10 -30 at k = to about 
10- 10 at k = 75. 

Let us turn now to the problem of restoring of the E(g) dependence at large g values. We shall 
focus on the ground state. First of all we note that in work || an interpolation expression for the 
ground (tightly bound) state energy was obtained. In spite of the multiple anticrossings at TC < 300 
and of the related computations complicating, the fit of || provides precision within 10 -3 -j- 10 -2 
in the range of TC values 0.1 < TC < 10 4 . 

The ground level energy asymptotic at large g (or, the same, at large TC) is given by 

E(H)^ 1 -H- 1 -\n'(cH) + ... (9) 

(see, for example, Here c is a dimensionless constant. First, we consider the possibility of 

restoring of the leading term parameters in (9) - the power index and the constant multiplier - 
using perturbation theory. Methods applicable to this problem are considered in || ||. Note that 
for coming of the asymptotic into action the leading term in (9) should be large comparing with the 
correction term. One can look, for example, at the results of work Jl3j (where the values of £(TC) 
were obtained by variational procedure) and ensure that only if TC > 10 2 then the binding energy 
will make less than 20% of So, we can speculate about restoring of asymptotic parameters 
only if we succeeded in summation of E(g) in this region of external fields. But we failed to do 
this having used only 75 coefficients Ek, so the linear asymptotic couldn't be restored. This was 
confirmed in our attempts to apply methods suggested in the works |J - no plausible result was 
obtained. In the method of Ref. ||, parameters of asymptotic of the function E(g) were linked to 
behaviour of coefficients Dk in dependence of their number k at large k. Namely, if E(g) — > Cg v 
at g — > 00, then in our case we get similarly to || 

Ck v ~ x 

Dfe ^<r(2z,)r(2z, + 6 )' 
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Then it was suggested to perform the fit of C and v using known Dk and their errors by means 
of the x 2 method. But in our case the value of x 2 i n hs minimum was extremely large (about 10 8 
even if we tried to fit only 5 coefficients Dk at statistical error o = 10 -10 , and we had no reason 
to increase this value of a). This result indicates that the asymptotic of Dk comes into action at 
values of k much larger than 75. 

The power index in E(g) asymptotic could be traced also using the method of ||. This method 

concerns of taking of the limit of the expression — — at w — > oo (or, the same, the limit of 

B(w) 

- 7^7 — — at y — > 1), that gives exactly the value of v. But numerical calculation showed that 

B i.y) 

in the region where B(y) was recovered (at y close to 1 we obviously should have increasing of error 
due to finite number of Dk used) we did not get reasonable precision for the limit value. 

Fig. 4 illustrates precision of the linear asymptotic. The curve plotted represents the binding 
energy (we used here formula (6) of Ref. 0) divided by Landau energy versus external field. One 
can see that at 7i < 100 error of asymptotic is more than 20%, and only at TC ~ 1000 precision 
reaches level of 1%. So, it appears to be impossible to obtain asymptotic parameters corresponding 
to Landau level for the Zeeman effect. 

Now, the question arises, whether we can subtract Landau energy from E(TC) and trace the 
second term of asymptotic. But this term doesn't work even at TC ~ 10 5 , which the graph of Fig. 5 
is to explain. We plotted there the values of binding energy, and the value of its logarithmic asymp- 
totic. We have chosen constant c in such a way, that the value of logarithmic term (1/2) ln 2 (c7^) 



of the asymptotic coincides with the data of [131 at TC — 10 . This occurs at c = 0.010 and the 



smallness of this constant indicates that the value of Tt is too small for speaking about asymptotic. 

And one can see that asymptotic curve and curve of exact data have considerably different 
slopes. It means that we have to recover dependence of the energy at the values of external fields 
larger than 10 5 . But it is really impossible having known only 75 coefficients E^. Thus, knowledge 
of 75 hypersucceptibilities did not allow to restore neither parameters of the Landau asymptotic 
nor logarithmic asymptotic of binding energy. 
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Figure 1: Binding energy for Is state in atomic units. The data evaluated in double precision (with 
help of Pade approximant [30/30] - solid curve , by straightforward summation - dotted curve) and 
in quadruple precision (with help of Pade approximant [30/30] - dashed curve). Crosses denote 
the data from Ref. @) 
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Figure 4: Relative precision of the linear asymptotic 
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Figure 5: Binding energy and its asymptotic. Solid curve is plotted using Eq. (6) from Ref. [BJ. 
Crosses denote the data from Ref. fl3 |. Dotted curve is the logarithmic asymptotic with c = 0.010 
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